arXiv: 1505.02028v3 [cond-mat.quant-gas] 27 Nov 2015 
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The absence of energy dissipation leads to an intriguing out-of-equilibrium dynamics for ultra¬ 
cold polar gases in optical lattices, characterized by the formation of dynamically-bound on-site and 
inter-site clusters of two or more particles, and by an effective blockade repulsion. These effects 
combined with the controlled preparation of initial states available in cold gases experiments can be 
employed to create interesting out-of-equilibrium states. These include quasi-equilibrated effectively 
repulsive ID gases for attractive dipolar interactions and dynamically-bound crystals. Furthermore, 
non-equilibrium polar lattice gases can offer a promising scenario for the study of quasi-many-body 
localization in the absence of quenched disorder. This fascinating out-of-equilibrium dynamics for 
ultra-cold polar gases in optical lattices may be accessible in on-going experiments. 


PACS numbers: 37.10.Jk, 67.85.-d, 03.75.Kk, 05.30.Jp 

Out-of-equilibrium dynamics of isolated quantum sys¬ 
tems has recently attracted a major interest H ig, in 
particular in the context of ultra-cold gases, where dissi¬ 
pation is basically absent [3]. Non-equilibrium quantum 
dynamics constitutes an exciting new field, notably in 
what concerns many-body localization (MBL), i.e. local¬ 
ization in excited states of interacting many-body sys¬ 
tems [4]. Recent cold-gases experiments are starting to 
unveil the non-trivial physics of MBL [5] . 

Although MBL is typically discussed in the presence of 
disorder, localization may occur in absence of it, as first 
discussed for ^He diffusion in ^He crystals [i[7]. Beyond 
a critical concentration, immobile ^He clusters could lead 
to percolation for the remaining ^He atoms. Quasi-MBL 
and glassy dynamics without disorder are attracting a 
growing attention, and various mechanisms for localiza¬ 
tion and eventual delocalization have been discussed 

0181 [IMS]. 

Meanwhile, experiments on magnetic atoms [IMHI 
and polar molecules mHni are starting to reveal the 
fascinating physics of dipolar gases. These gases are 
markedly different from their non-dipolar counterparts 
due to the long-range anisotropic character of the dipole- 
dipole interaction (DDI) [22l |23] . Polar gases in optical 
lattices (OLs) offer exciting possibilities for the study of 
lattice models [23] and quantum magnetism [24l |25] . 

In this Letter, we study non-equilibrium dynamics 
of ID polar lattice gases. This dynamics is character¬ 
ized by dynamically-bound on-site and inter-site clus¬ 
ters (BCs) generalizing on-site repulsively-bound pairs in 
non-polar gases [261128] . and by blockade repulsion (BR). 
We show how these effects result in interesting out-of¬ 
equilibrium states, including repulsive gases with attrac¬ 
tive DDI and dynamically-bound crystals. Moreover, po¬ 
lar lattice gases allow for quasi-MBL without disorder, as 
we illustrate for the set-up of Fig.[^ These scenarios can 
be realized in current experiments on polar molecules in 
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FIG. 1: (Color online) (a) A dimerized lattice can be used to 
create a gas of singlons and dynamically-bound NN dimers; 
(b) due to the BR a NN dimer (singlon) forms a block D = 
0011 (S = 001); (c) effective lattice formed by the blocks D, S, 
and additional empty sites. This effective lattice is employed 
in model § to show the realization of quasi-MBL. 

OLs. 

Model- We consider polar bosons in a ID OL. For 
a deep lattice the system is described by the extended 
Bose-Hubbard model (EBHM) [22j [23] 

= + + r ^ (1) 

(ij) i i,r>0 

where (• • •) denotes nearest neighbor (NN), bi (6t) de¬ 
stroys (creates) bosons at the i-th site, rii = b\bi, J is the 
hopping rate, U characterizes the combined on-site short- 
range interactions and DDI, and V/r^ is the strength of 
the DDI between sites placed r > 1 sites apart. J, U and 
V can be tuned independently by changing the lattice 
depth, the transverse confinement [29], the orientation 
and strength of the polarizing field {V < 0 for polariza¬ 
tion along the lattice axis), and by Feshbach resonances. 

Bound pairs.- We revisit first the concept of bound 
pairs in non-polar gases (V = 0). Doubly-occupied sites 
are characterized by an interaction energy U. If \U\ ^ J, 




2 



K/n r 


FIG. 2: (Color online) (a) Energy spectrum as a function of 
the center-of-mass quasi-momentum K for two particles in 40 
sites for = 0 and V = — lOOJ; density of the bound-states 
as a function of the relative coordinate r for the three bound 
states at K = 0, as marked by the blue circles in Fig. (a), at 
energies E/J ^ —5 (b), —13 (c) and —100 (d). 

energy conservation maintains on-site pairs irrespective 
of the sign of U [30]. For U > 0 those pairs, also 
called repulsively-bound pairs m, are hence dynamically 
bound. Conversely, two separated particles cannot be 
brought to the same site, i.e. singlons experience hard¬ 
core repulsion m- However, singlons may resonantly 
move through on-site pairs since a single-particle hopping 
swaps doublon and singlon positions (21 ^ 12) [3^ [33] . 

The long-range DDI allows for dynamically-bound 
inter-site pairs [34|. Figure a) depicts a typical two- 
particle spectrum, for 1/ = 0 and V = — lOOJ. For each 
center-of-mass quasi-momentum K G [—7r,7r], the spec¬ 
trum presents a continuum of scattering states and a dis¬ 
crete set of isolated inter-site bound states (BSs) [351(36] . 
which as for on-site bound pairs in non-polar gases, are 
maintained by energy conservation, irrespective of the 
sign of V. Figures [^b-d) show the probability of finding 
two particles r sites apart for the BSs at K = 0. For 
binding energies close to the continuum, the relative po¬ 
sition of the pair delocalizes over many sites (Fig. |^b)). 
These delocalized BSs are for any practical purposes in¬ 
distinguishable from the scattering states. Instead, as 
shown in Fig. [^c-d), deeper BSs present a well defined 
relative distance r. Below we restrict the term bound 
pair (BP) to deep BSs at fixed r < Tc, where the criti¬ 
cal Tc is defined as the largest r satisfying the condition 
/(r) = 2{JIV)‘^l{r~^ — (r + 1)~^)^ 1 [37]. Note that 

even if U = 0, the inter-site DDI stabilizes an on-site BP 
that is buried within the scattering states in Fig. |^a), 
close to energies E = U = 0. On-site and NN BPs de¬ 
mand \V — U\ ^ J to avoid resonances between on-site 
and inter-site interactions [38] . 

Bound clusters.- On-site interactions may bind more 
than two particles in on-site BCs. However, on-site BCs 


are unstable against three-body losses and play a relevant 
role only at relatively large lattice fillings [39] . Polar lat¬ 
tice gases allow for inter-site BCs of more than two par¬ 
ticles, each particle being within a distance r < Vc of at 
least another particle of the cluster. Several important 
points must be noted. First, although sites with more 
than one particle may be involved, inter-site BCs are typ¬ 
ically formed by singly-occupied sites, and hence these 
clusters are in general stable against three-body losses. 
Second, whereas on-site BPs are obviously precluded for 
polarized Fermi gases, inter-site BPs and BCs are possi¬ 
ble even in that case. Third, in contrast to on-site BCs, 
inter-site BCs may present internal resonances, e.g. the 
cluster 1101 may remain bound, but resonates with 1011. 
BCs are a general feature of non-equilibrium polar lat¬ 
tice gases in any dimension even at low fillings, as long 
as the DDI is large enough. In particular, massive BCs 
of a size comparable to the whole system may be formed 
if the mean-interparticle distance R < Vc. An example 
of massive BC is provided by particles initially placed at 
regular distances Vin < Tc- The absence of dissipation 
maintains this dynamically-bound crystal (Fi gia)). 

Blockade repulsion.- The formation of inter-site BPs 
has as a counterpart a vanishing probability of finding 
the particles at a distance r < Tc in loose inter-site BSs 
and scattering states (rc = 2 for Fig. [^b) [37]). This 
exclusion region leads to an effective BR between parti¬ 
cles initially at a distance > Tc- This BR becomes 
evident in the density-density correlation g 2 {t,r) = 

{ni{t)ni^r{t))/{{ni{t)){ni^r{t))). If Gn > Tc at t = 0, 
the subsequent dynamics shows BR, i.e. g 2 {t > 0,0 < 
^ < ^c) = 0, as discussed below. 

Repulsive gas for attractive DDL- Combining BR 
with a proper initial-state preparation allows for the cre¬ 
ation of a repulsive gas for attractive DDI. Such a gas 
may be realized by placing particles at the minima of 
a superlattice with period Vin > Vc and subsequently 
removing the superlattice; atoms in sites with more 
than two particles may be eliminated by using resonant 
light [40]. Under these conditions no BP or BC is present, 
and the system forms a singlon gas with effective BR at 
radius Tc- We have performed time-dependent density- 
matrix renormalization group (t-DMRG) simulations [41] 
for V = —lOOJ and U = 0 {vc = 2) and 4 particles ini¬ 
tially Vin = 5 sites apart. After a short time ^ J“^, 
the density (uj) and g 2 {r = 0) converge to the values 
expected for a homogeneous gas (inset of Fig.|^b)) [42] . 
However, due to BR, ^ 2 (CO < ^ '^c) = 0 for all t, 

whereas ^2(C ^ > ^c) bas a non-trivial dynamics reaching 
quasi-equilibrium (Fig. |^b)). The time-averaged single¬ 
particle correlation gi{r) = ill shown 

for different times ti ^2 in Fig. [^c), also indicates quasi¬ 
equilibrium. The equilibration of gi and ^2 relies on the 
absence of BPs or BCs, contrasting with the quasi-MBL 
scenario below. 
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FIG. 3: (Color online) (a) Dynamically-bound crystal: density (ui) (obtained using t-DMRG) for U = 0 and V = —lOOJ (rc = 
2), for 4 particles in 24 sites initially placed Vin = 3 sites apart. Note that since Vin = rc + 1, there is still some residual 
dynamics. For smaller Tin < rc the crystal is perfectly preserved in this time scale, (b)-(c) Effective repulsive ID gas: same 
parameters as in Fig. (a) but Vin — 5: (b) g 2 {t,r > 0) at different times (in the inset we show g 2 (t,r = 0) |l2]) and (c) ^i(r), 
for J{t 2 — ti) = 5 and different t 2 values. For both gi and ^ 2 , r denotes the distance from an initially occupied site. Note 
that not only the on-site density and density fluctuations converge to quasi-equilibrium, but also gi and g 2 at longer distances. 
Note as well that g2{t,0 < r < Vc) = 0 due to BR. 


Although the effective repulsive ID gas resembles a 
super-Tonks gas [43], the physics behind is very differ¬ 
ent. In the super-Tonks case, an initially repulsive gas 
is dynamically brought into an attractive regime. Even 
if in that regime the two-body ground-state is a bound 
state, in absence of dissipation the system remains in an 
excited state characterized by inter-particle repulsion. In 
contrast, BR is crucially maintained by both the absence 
of dissipation and by the lattice, which provides a finite 
band- width and discrete particle motion. 

Quasi-many-body localization.- Polar lattice gases of¬ 
fer interesting possibilities for the study of quasi-MBL 
without disorder. BCs of M particles move as a whole 
with hopping and hence BCs with M > 1 

are for any practical purposes immobile (although in¬ 
cluster quasi-resonances may be still possible). As for 
^He (6] [7], massive BCs and BR may induce percola¬ 
tion for large-enough filling and \V\/J. Interestingly, as 
shown below, ID polar lattice gases may present quasi- 
MBL even for low fillings {R ^ Vc) and moderate DDI 
achievable in experiments. 

We illustrate the possibilities of ID polar lattice gases 
for quasi-MBL within a simplified scenario. Resem¬ 
bling the recent experiment of Ref. [5] we consider a 
dimerized OL such that only the lower sites are pop¬ 
ulated (Fig. [^a)). After eliminating atoms in doubly- 
occupied sites [40], neighboring lower sites are hence ei¬ 
ther both occupied, or only one of them, or none. Then 
the superlattice is removed. We consider \V\/J such that 
Tc = 2, and hence NN dimers form a BP, and a BR at 
Tc = 2 is established. As a result, blocks 0011 and 001 
behave as well-defined particles that we call D and S, re¬ 
spectively (Fig. [^b)). We neglect DDI for r > Vc neigh¬ 
bors, since it is well within the bandwidth, and obtain 


the effective model [44]: 

Eff = - E (A^^i + JdDIE + ^Dls]SiDfj , (2) 

i'i'j) 

where i, j denote the sites of the effective lattice formed 
by D’s, S’s, and empty sites of the original lattice belong¬ 
ing neither to a D or an S (Fig. s c)). In H,ff, Dj (Sj) 
destroys a D (S) at the effective site j. Assuming for 
simplicity U ^ V, J, the D hopping rate is Jd = 

The third term in Heff is the swap DS ^ SD, occur¬ 
ring at rate Note that a site of the effective 

lattice is occupied by a D, an S, or empty. Due to this 
hard-core constraint, in absence of swapping, D’s and S’s 
would trivially localize each other, for any ratio Jd!J- 

Swapping allows for the motion of D’s and S’s. How¬ 
ever, the fact that ^ J for J/\^\ ^ 1 rnay result 

in localization following similar arguments as in Ref. [3]. 
For J/\V\ 0, the motion of S’s is blocked by D’s. For 

finite J/\V\ <C 1, the motion of D’s changes the energy of 
the S gas [45]. If this change, AE ^ the motion of 

D’s is hindered. However, limited quasi-resonant D mo¬ 
bility, involving AE < remains possible, leading 

to partial D diffusion at times ^ l/H. 

We have performed exact diagonalization calculations 
with periodic boundary conditions (PBC) of the evolu¬ 
tion of the many-body state |T(t)) given by Model (© 
for small systems {ud^tls^L) of ud D’s, Ug S’s and L 
lattice sites, corresponding to Le// = L — Sud — 2ns 
effective sites [44]. We average over various initial 
random distributions of D’s and S’s at fixed positions 
in the effective lattice. Figure shows for = 

0.013 iy = — lOOJ) the dynamics of the inhomogene- 
ity of D’s, AAT ^ - y+i|^(f))P 

(with y ^ f)]t)^) [3]. Perfect homogeneity means 
A A = 0. In Fig. [^ we depict for comparison the re¬ 
sults for Ei/J = 0.13 [46]. Whereas for Ei/J = 0.13, D’s 
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FIG. 4: (Color online) Inhomogeneity AN{t) as a function 
of Qt for Q/J = 0.013, for {nD,ns,L) = (2,2,20) (dashed 
blue) and (3,3,27) (solid red). The results are obtained by 
exact diagonalization averaging over 50 and 25 random ini¬ 
tial conditions, respectively. The shaded region indicates ap¬ 
proximately the region of fast decay due to quasi-resonances. 
We depict for comparison the results for (3,3,27) for Q/J = 
0.13 (dotted pink). The inset shows the IPR for (3,3, 27) for 
Q/J = 0.013 (solid red) and Q/J = 0.13 (dotted pink). 

diffuse within a time scale l/f^, for Q/J = 0.013 quasi¬ 
resonances allow only a fraction of D’s to delocalize in 
this time scale (shaded region in Fig. and a much 
slower dynamics follows. This slow dynamics is charac¬ 
teristic of systems with PBC due the collective motion 
of all D’s [4]. Consistent with this, the time scale of the 
slow D dynamics for the {ud^tis^L) = (3,3,27) case is 
approximately 10 times longer than that for ( 2 , 2 , 20 ) [47] . 
We hence expect an exponentially diverging time scale for 
the slow-dynamics for growing number of dimers. 

We may expand |^(t)) = ^)l^) 

nmax = many-body states |z/) account- 

ing for all possible distributions of S’s and D’s in the 
effective lattice. MBL may be visualized as localiza¬ 
tion in this many-body space. The latter is best quan¬ 
tified by the inverse participation ratio (IPR), r]{t) = 

^ fully delocalized (localized) 
state presents 77 ^ 1 (^ l/'^macc)- Whereas for 0/J = 
0.13, r]{t) ^ 1 at Dt < 10, for VtjJ = 0.013, r](f) remains 
very small even for Qt ^ 1 (inset of Fig.[^, showing the 
appearance of quasi-localization in the many-body space. 

Experimental feasibility. - The previous scenarios can 
be realized with polar molecules in OLs, as we illus¬ 
trate for the case of NaK, which possesses an elec¬ 
tric dipole of 2.72 Debye in its lowest ro-vibrational 
level [48]. We consider the realistic case of partially 
polarized molecules with 1 Debye. For a lattice spac¬ 
ing of 532nm, V/h IkHz. Assuming a lattice depth 
of l^Ereci with Er^f^lh 2.75 kHz the recoil energy, 
J/h 10 Hz= IH1/100, and hence 1 /D ^ Is. As shown 
in Fig. [^the slow dimer dynamics can be much larger. 


stretching well beyond a minute, which is the typical 
maximal life time in experiments. Localization can be 
explored either in expansion experiments, or by site- 
resolved measurements. Moreover, the formation of a 
repulsive gas with attractive DDI can be readily moni¬ 
tored. Tightening an overall harmonic trap should result 
in the formation of an incompressible crystalline core that 
can be revealed by measuring the saturation of the mean 
radius of the sample and/or by site-resolved measure¬ 
ments. Furthermore, BR hinders two or more molecules 
to gather at the same site, preventing chemical recombi¬ 
nation losses despite attractive DDI. 

Summary.- The absence of dissipation leads to rich 
out-of-equilibrium dynamics in polar lattice gases char¬ 
acterized by the formation of inter-site bound clusters 
and blockade repulsion even for attractive DDI. The 
combination of these effects with the control possibili¬ 
ties of ultra-cold gases may allow the realization of effec¬ 
tive repulsive ID gases with attractive DDI, the creation 
of dynamically-bound crystals, and most interestingly, 
quasi-MBL in absence of disorder. The latter opens in¬ 
teresting perspectives for observing a dynamical phase 
transition in polar lattice gases from a delocalized to a 
quasi-MBL regime as a function of the V/ J ratio. 
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0.13 does not describe a polar gas at Y = — lOJ, which 
has rc = 1, and hence cannot be mapped to model §. 

[47] Gollective motion of two D’s occurs in a time scale 
T 2 ~ AE/D^, where AE is the typical energy shift when 
displacing one D. In our simulations Q.T 2 — 10. Eor the 
case of three D’s the time scale for the collective motion 
is T 3 ~ AE^/D^, and hence Qt 3 ~ 100 in good agreement 
with our results. 

[48] J. W. Park, S. A. Will, and M. W. Zwierlein, 
arXiv: 1505.00473 
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SUPPLEMENTARY MATERIAL FOR 
’’OUT-OF-EQUILIBRIUM STATES AND 
QUASI-MANY-BODY LOCALIZATION IN 
POLAR LATTICE GASES” 

We consider atoms occupying sparsely the lower sites 
of the lattice of Fig. 1(a) of the main text. Neighboring 
lowest sites may be both occupied, one occupied and one 
empty, or both empty. Atoms in doubly-occupied sites 
or in sites with even higher occupation are removed by 
resonant light. As a result, when the additional superlat¬ 
tice is turned off, one has either nearest-neighbor (NN) 
dimers, or singlons (i.e. singly-occupied sites) that are 
at least three sites apart from the nearest singlon or NN 
dimer. This provides our initial condition at t = 0. 

We assume a ratio V j J such that rc = 2^ ensuring that 
this minimal distance persists during dynamics. Namely, 
due to blockade repulsion no particle can ever be found 
at a distance equal or less than two sites from a singlon 
or a NN dimer. As a result each singlon and each doblon 
have always two empty sites at their left and right. Hence 
we can safely define the blocks S' = 001 for singlons and 
D = 0011 for NN dimers [1]. Due to periodic boundary 
conditions, any possible state maps into a distribution of 
S’s and D’s, and additional empty sites 0, which do not 
belong to any S or D. Therefore, we map the real system 
onto an effective lattice with sites that are either empty, 
occupied by an S, or occupied by a D (see Figj^a)). 

The tunneling of S’s to empty effective sites takes place 
via single-particle hopping J . D’s can also move of one 
effective site via second-order processes in two ways: (i) 
oil ^ 020 ^ 110, with an amplitude J‘^/{U — V); 
and (ii) 011 ^ 101 ^ 110 (see Fig. mb)) with an 
amplitude SJ‘^ jlV. Hence the doublon-hopping rate is 
Jd = /{U — V) SJ‘^f7V. Assuming, for simplicity, 

U ^ y, J, we approximate Jd = SJ^/TU, as in the main 
text. Note however that this assumption is not strictly 
needed. It is however necessary, as we mention in the 
main text, to fullfil the condition \U — V\ ^ J. Other¬ 
wise NN dimers become mobile at a hopping rate similar 
to J (exactly equal to J in the case U = U) [2]. Fi¬ 
nally, neighboring D’s and 5'’s may swap their position 
via the second-order process sketched in Fig. [^c), with 
an amplitude Q = 

For Tc = 2 and due to the 1/r^ decay of the dipole- 
dipole potential, the interaction at distances larger than 
two sites is, by definition of Tc, well within the band¬ 
width, and can hence be neglected in a good approxima¬ 
tion. As a result, the dynamics of the polar lattice gas 
reduces to the hopping of 5'’s and D’s, and the swap of 
D’s and 5'’s; this dynamics is described by model (2) of 
the main text. 

Note that this effective model resembles the one previ¬ 
ously introduced for binary mixtures of light and heavy 


particles ISHS] in the context of MBL. However in [SHS, 
heavy particles block light particles, but light particles 



FIG. 5: (Color online) (a) The outer ring represents a ID 
polar lattice gas with periodic boundary conditions, under 
the conditions discussed in the text. Singlons form effective 
clusters S' = 001 (orange) and NN dimers effective clusters 
D = 0011 (green) [T]. The inner ring represents the effective 
lattice corresponding to the real lattice in the outer ring . 
Each site of the effective lattice may be either empty (white), 
occupied by an S (orange), or by a D (green). Figures (b) and 
(c) depict, respectively, the second-order processes responsible 
for the hopping of D’s into a neighboring empty effective site, 
and for S-D swaps (see text). 


do not block heavy ones. In model (2), D’s block S’s 
and vice versa, but S’s and D’s can swap their positions. 
Note that, as mentioned in the main text, these swaps are 
crucial, since without S-D swaps (and as a consequence 
of the hard-core constraint of the effective model), D’s 
and S’s would block each other mutually, and the system 
would be trivially localized. Indeed, it is the fact that 
the particles can in principle extend over the whole lat¬ 
tice via S-D swaps that makes the problem non-trivial. 
As a consequence, the crucial parameter that controls the 
quasi-many-body localization is Q/J, as shown in Fig. 4 
of the main text. 


[1] Obviously we could have chosen D = 1100 and S = 100, 
and the mapping to the effective model would be identical. 

[2] Resonances between on-site interaction U and NN interac¬ 
tion V may lead to interesting physics, which lies beyond 
the scope of the present analysis. 
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